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Abstract 

The degenerate parabolic equation u t + d x [u 3 (u xxx + u x — sinx)] = models the evolution 
of a thin liquid film on a stationary horizontal cylinder. It is shown here that for each mass there 
is a unique steady state, given by a droplet hanging from the bottom of the cylinder that meets 
the dry region with zero contact angle. The droplet minimizes the associated energy functional 
and attracts all strong solutions that satisfy certain energy and entropy inequalities, including all 
positive solutions. The distance of solutions from the steady state cannot decay faster than a power 
law. 

Keywords Thin liquid film; coating flow; strong solutions; steady state; symmetrization; energy; 
entropy method; Lyapunov stability; power-law decay. 

AMS Subject Classification: 35K25; 35K35; 35Q35; 37L05; 76A20. 

1 Introduction and description of the results 

Degenerate fourth order parabolic equations of the form 

Ut + V • (u n VAu) + lower order terms = (1.1) 

are used to model the evolution of thin liquid films on solid surfaces. Here, u(x,t) describes the 
thickness of the fluid at time t at the point x, the fourth derivative term models the surface tension, 
and the exponent n > is determined by the boundary condition between the liquid and the solid. 
The equations were derived from the underlying free-boundary-value problem for the Navier-Stokes 
equation by the lubrication approximation, which is valid if the film is relatively thin. 
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Figure 1 : Thin liquid film on the outer surface of a horizontal cylinder. 



An interesting example is Pukhnachev's model for a thin liquid film on a rotating horizontal cylin- 
der 11331 . as shown in Figure 1. The model equation is 



u t + d e 
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Here, 9 e K/ (2ttZ) denotes the angle measured from the bottom of the cylinder, the film is assumed to 
be uniform in the axial direction, and inertia is ignored. The power n = 3 indicates a no-slip boundary 
condition between the liquid and the solid. The other physical parameters of the system are the surface 
tension (a), the viscosity (u), the acceleration of gravity (g), the density of the fluid (p), the radius of 
the cylinder (R), and the rotation speed (a;). Note that the surface tension appears with an additional 
lower-order term that corrects for the curvature of the cylinder. Eq. (11.21) refines an earlier model 
of Moffatt ll30l that neglects surface tension. Numerical and asymptotical analysis of Pukhnachev's 
equation along with numerous open questions can be found in [f2l |3l |2TJ , and linearizations about 
steady states are examined analytically and numerically in |[T3~lfT9l . 

We will show that the long-time behaviour of Pukhnachev's model on a non-rotating cylinder is 
controlled by steady states. Our results imply that 

• Eq. (11.21) with to = has for every given mass a unique nonnegative steady state; 

• the steady state minimizes energy and attracts all solutions of finite entropy; 
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• the distance of solutions from the steady state decays no faster than a power law ~ t~s. 

Note that solutions of finite entropy are almost everywhere strictly positive. The steady state has the 
shape of a shallow drop hanging from the bottom of the cylinder, with a dry region at the top that it 
meets at zero contact angle, see Figure 2 below. We suspect that the distance from the steady state 
actually behaves as t~3. This conjecture is supported by simulations, and by analogy with aggregation 
processes such as late-stage grain growth in alloys 11271 and the formation of drops at a faucet, where 
the growth of the grain or droplet is limited by the rate of mass transfer through a region of low density 
to the region of accumulation. 

Our results are motivated by the work of Carrillo and Toscani, who proved global convergence to 
self-similar solutions for the thin film equation 

u t + d x [u n u xxx ] = (1.3) 
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on the real line with n = 1 [[Toll . In contrast with the recent precise convergence results of Giacomelli, 
Kniipfer and Otto lfT6l . our conclusions are qualitative and global. We work in a large class of non- 
negative strong solutions that contains all positive classical solutions, and we do not assume (and do 
not prove) uniqueness of the solutions. Our lower bounds on the distance from the steady state should 
be compared with results of Carlen and Ulusoy |fTT|. who showed for Eq. (|1.3I) with n — 1 that the 
distance from the self-similar solution satisfies a power-law upper bound. We will give a more detailed 
description at the end of this section. 

Thin liquid films have been the subject of rigorous mathematical analysis since the pioneering 
article of Bernis and Friedman [5 J. A vast body of papers is dedicated to existence of solutions, 
regularity, long-time behaviour, finite-time blow-up, and the interface between wet (u > 0) and dry 
(u = 0) regions, see for example (TJ H U\ HI CHI [3TJ and references therein. An even larger part of 
the literature studies the properties of physically relevant solutions through asymptotic expansions, 
numerical analysis, and laboratory experiments. 

A fundamental question is well-posedness: which initial values give rise to unique nonnegative 
solutions that depend continuously on the data? The difficulty is that solutions of fourth order parabolic 
equations generally do not satisfy a maximum principle, and linearization leads to semigroups that do 
not preserve positivity. To give a simple example, the function u(x, t) — 1 + 1 cos(x), which develops 
negative values after t — 1, solves Eq. (11.41) (given below) with n = and a = 1. Bernis and Friedman 
proved that initial-value problems in one space dimension have weak solutions in suitable function 
spaces 0. A far-reaching technical contribution was their use of energy and entropy functional that 
decrease along solutions. Still, after twenty years, uniqueness remains an open problem. 

There are many questions surrounding steady states: Are they uniquely determined by their mass, 
and if not, how many are there? Are they strictly positive, and if not, what is the contact angle 
between wet and dry regions [|25l ? Under what conditions are they stable, do they attract all bounded 
solutions, and what is the rate of convergence? Since energy decreases along solutions, we expect that 
steady states should correspond to critical points of the energy, that solutions should converge to steady 
states, and that minimizers of the energy should be asymptotically stable. However, in the absence of a 
proper well-posedness theory, the proof of these statements requires more than a standard application 
of Lyapunov's principle (as stated, for example, in lfT8l ). 

One strategy for proving convergence to equilibrium is to use entropy in place of energy. The 
basic thin film equation (11.31 ) has a family of entropy functionals of the form Sp(u) = J n dx that 
decrease along solutions. For (3 = n — 2, this was established by Bernis and Friedman, for (3 = n — | 
it is due to Kadanoff (see [6j), and the range n — 3 < (3 < n — | was developed independently in |@]|7]|. 
Other families of entropies have since been discovered EOl 1241 . In their classical papers, Beretta, 
Bertsch, Dal Passo and Bertozzi, Pugh used these entropies to show that solutions of Eq. (11.31) with 
n > on an interval become positive after a finite time and converge uniformly to their mean [@]|7]|. 
Remarkably, this convergence holds for a very broad class of weak solutions, about which little else 
is known. Several works over the last decade have combined energy and entropy methods by deriving 
coupled inequalities for the energy and entropy dissipation. In this way, Tudorascu proved exponential 
convergence to the mean for thin films on finite intervals 11351 . and Carrillo, Toscani [10J and Carlen, 
Ulusoy ifTTTl have proved global (power-law) convergence to self-similar solutions of the thin equation 
with n = 1 on the real line [ITOl ITT|| . 

A promising approach, first pursued by Otto [f3~T| , treats thin-film equations as gradient flows that 
evolve by steepest descent on the space of nonnegative functions of a given mass, endowed with a 
suitable metric. Since Eq. (11.11) has the form u t = [u n (§f jj , a formal computation suggests certain 
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Figure 2: The energy minimizing steady state for Eq. d 1 -41) with a = 1 and initial data u = 0.5, 1, 2. The 
minimizer is concentrated if the mass is small, and spreads out over (— n, ir) as the mass increases. 



variants of the Wasserstein metric. These metrics measure the distance between two mass densities 
as the cost of an optimal transportation plan that pushes one density forward to the other. Gradient 
flow methods have been very successful for the porous medium and fast diffusion equations, which 
are second-order degenerate parabolic equations that share many features of thin films. There also has 
been some progress on thin film equations with n = 1, where the metric is just the standard optimal 
transportation distance OT1 . Thin-film equations with n > 1 have proved resistant to this approach, 
because the geometry of the relevant Wasserstein distance with mobility is not well understood, and 
the thin film energy is not geodesically convex 0. 

Yet another strategy is to linearize about the steady state. This linearization is delicate if the steady 
state has a dry region. Recently, Giacomelli, Kniipfer and Otto have developed a technique specifically 
for droplets with zero contact angles ifToll . They prove well-posedness and convergence to the steady 
state for initial values in a singularly weighted Sobolev space that forces the solution to vanish to a 
certain order at the contact point. So far, these techniques have been developed for Eq. (11.31) with 
n = 1. 

Outline of the paper. We study the role of energy-minimizing steady states for the dynamics of 

u t + d x [u n [u xxx + a 2 u x - sins)] = , igS1 = R/(2ttZ) (1.4) 

with n > and a > 0. Each of these equations describes the evolution of a thin film. As in Eqs. (11.11) 
and (|1.3I) . the exponent n is related to the boundary condition between the liquid and the solid, the 
first two terms in the parentheses model surface tension, and the third term accounts for gravitational 
drainage. The o?d x [u n d x u] term is reminiscent of the porous-medium equation, but appears with the 
opposite sign, resulting in a long-wave instability 10. The units of time and length are scaled so 
that the surface tension and gravitational terms appear with coefficient one. The parameter a > is 
a geometric constant, with a = 1 for a horizontal cylinder. Values of a > 1 appear when inertial 
effects are taken into account [|23l Eq. (2.3)]. In general, the geometric coefficient is a function that 
depends on the curvature of the surface [29, Eqs. (64) and (68)]. We are mostly interested in the case 
where n = 3 and a = 1, which corresponds to Pukhnachev's model at zero rotation speed, but find it 
illuminating to consider also other values of n and a. 
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Many competing definitions of weak solutions have been proposed. Following Bernis and Fried- 
man, we define a strong solution of Eq. (11.41) on a finite time interval (0, T) to be a nonnegative function 
uGL 2 ((0,T), H 2 (n)) that satisfies 

J J lyU(j) t — (u-xx + c? u + cosx)(-u n 1 .) a .| dxdt = (1.5) 

for every smooth test function (ft with compact support in (0, T) x Vt. Such solutions are believed to 
be unique. We will consider strong solutions that exist for all t > and satisfy additional bounds on 
the energy and entropy (see Section B}. In particular, these solutions are almost everywhere strictly 
positive. 

Our main results concern the convergence of solutions to steady states. We start from the energy 
for Eq. (11.41) . given by 

E(u) = - I (u 2 — a 2 u 2 ) dx — / ucosxdx, (1.6) 
n Jn 

where the first integral accounts for the surface tension, and the second integral for the gravitational 

dt 
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he surface tension, and 
potential energy. Formally, the energy decreases according to dE j^ = —D{u), where 



D(u) — / u n (u xxx + a 2 u x — sin x) 2 dx (1.7) 

J{u>0} 

is the dissipation associated with Eq. (11.41) . 

Note that the energy is not convex for a > 1, and not bounded below in H 1 for a > 1. Nev- 
ertheless, we show in Section [2] that for every choice of a > the energy has a unique minimizer 
among nonnegative ^functions of a given mass M. We denote this minimizer by u*. The shape of 
u* (depending on a and the mass) is described precisely in Theorem 12 .4[ The minimization problem 
is complicated by the fact that the minimizers can have dry regions, where the Euler-Lagrange equa- 
tion is weakened to an inequality. We combine a careful analysis of the variational inequality with 
symmetrization techniques. 

In Section [3l we show that u* is the unique steady state of Eq. (11.41) with zero dissipation when 
a < 1. For a > 1, we find also saddle points, as well as a continuum of two-droplet steady states that 
are not critical points for the energy in L 2 . Since the definition of strong solutions forces steady states 
to meet any dry spot at a zero contact angle, Theorem 12.41 implies in particular that such steady states 
exist if M(l — a 2 ) < 2n. For each value of a and M, there is also a continuum of time-independent 
solutions of Eq. (11.41) with compact support and positive contact angles, which are analogous to the 
steady states in E51 . Their role in the long-term evolution of positive solutions remains open. 

Section ID contains the main convergence result, Theorem 14. 31 We show that the energy minimizer 
u* is a dynamically stable, locally attractive steady state of Eq. (11.41) . For a < 1, we show that 
all solutions that satisfy certain energy and entropy inequalities converge to u*. For a > 1, the 
convergence holds on a sub-level set of the energy that contains no other strong steady states. 

The entropy methods of flU U\ SSI [HI ESI do not apply to Eq. (11.41) . because the entropy can 
increase as well as decrease along solutions. For steady states with dry regions, the entropy is not even 
finite. We combine Lyapunov's method with a linear bound on the growth of the entropy to produce a 
sequence of times along which the solution converges weakly to a steady state, using a recent argument 
of lfT2l . We pass to convergence in norm along the full solution by proving a local coercivity estimate 
on the energy near the minimizer (which becomes global for a < 1). 
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Finally, in Section [51 we turn to the rate of convergence. We show that convergence to the steady 
state is exponential whenever u* is strictly positive. This extends the results of 0] |7j [35J to examples 
where the steady state is non-constant (but note that our proof uses energy in place of entropy). When 
u* has a dry region and n > |, we show that solutions cannot approach it faster than a power law 

2 

~ t 2 ™- 3 . In the proof, we show that Kadanoff 's entropy S n _3 grows at most linearly. 

All our results are easily adapted to the long-wave stable case of Eq. ( 11.41) where the a 2 d x [u n d x u] 
term appears with the opposite sign (see [|8l). There, the energy-minimizing steady state is dynamically 
stable and attracts all solutions of finite entropy. It is strictly positive and exponentially attractive so 
long as M(l + a 2 ) > 2ti. Otherwise, the rate of convergence is limited by a power-law, at least when 
n > |. For M (1 + a 2 ) = 2tt, the minimizer has a touchdown zero, and for M(l + a 2 ) < 2ir, it has 
the shape of a droplet with zero contact angles. 

Acknowledgments. This work was partially supported by NSERC though a Discovery grant 
(A.B., B.S.) and a postdoctoral fellowship (M.C.). We thank Victor Ivrii, Mary Pugh, and Roman 
Taranets for stimulating discussions, and Rick Laugesen and Vladislav Vasil'evich Pukhnachev for 
valuable comments on an early version of the manuscript. 



2 Energy minimizers 

In this section, we study the energy landscape over the space of nonnegative functions of a given mass, 

C M = |m g L 2 {Vt) | u > 0, J u{x) dx = ikfj . 

We use two topologies on this space, the usual L 2 -distance \\u — v\\ 2 , and the -ff 1 -topology, with 
distance function 

d H i(u,v) = \\u x - v x \\ 2 , u,veC M - 

Since u and v have the same mass, their difference has mean zero, and this distance is equivalent to 
the usual i/ 1 -distance. It is clear from Eq. (11.61) that the energy is continuous in H 1 . By convention, 

E(u) = ooifw £ H l . 

For a < 1, E is a positive definite quadratic form, and hence strictly convex. This can be seen 
either from the Wirtinger inequality (see, for example, [fT5l p. 61]), or by writing the energy in terms 
of the Fourier series of u as 

M 2 

E(u) = 7T J2(P 2 - * 2 )WP)\ 2 - o?— - tt(«(1) + «(-l)) . 

p&\0 

For a = 1, the energy is convex, but not strictly convex on Cm, because its Fourier expansion depends 
linearly on w(±l). For a > 1, convexity is lost. 

We will show that the energy has a unique minimizer on Cm, and describe its profile. Our first 
lemma shows that a minimizer exists. 

Lemma 2.1 (Existence of minimizers.) For every M < oo, E attains its minimum on Cm- 
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Proof. Let u E Cm- Since u is nonnegative and has mean it satisfies 

/ u 2 dx < M \\u\\ L oo , \\u\\ L oo < — + y/ir\\u x \\ L 2 . 
Jn ^ 

Inserting these estimates into the functional, we obtain 

E(u) > - (\\uzWi? £-) - — W - (l + -)M , 

which shows that E is bounded below on Cm- 

Consider a minimizing sequence {%}j>i- By Eq. (12.11) . the sequence is bounded in H 1 . We invoke 
the Rellich lemma and pass to a subsequence (again denoted by {«,}) that converges weakly in H 1 
and strongly in L 2 to some function u* in Cm- Since E is weakly lower semicontinuous on H 1 , we 
have 

inf E(u) < E(u*) < lim E(uj) < inf E{u) , 
and conclude that E attains its minimum at u*. □ 
The Euler-Lagrange equation for E under the mass constraint is given by 

u xx + c?u + cos X = X , (2.1) 

where A is a Lagrange multiplier. We need to incorporate also the positivity constraint. If u 6 Cm, we 
decompose Vt according to the value of u into the positivity set and the zero set of u, defined by 

p(u) = {x e Q | u{x) > 0} , Z(u) = {x e Q | u(x) = 0} . 

Lemma 2.2 (Euler-Lagrange equation and zero contact angle.) If u* minimizes E on C M , then it 
solves ( 12. i| ) on P{u*). The Lagrange multiplier is positive and satisfies A > supjcosa; | x e Z(u*)}. 
Furthermore, u* is of class C 1,1 , andu* x = on dP{u*). 

Proof. Let be a smooth 2n -periodic test function, and set 

M f 

u £ = jy-{u* + e<p) , M £ = (u* + e<f))dx. 



We compute the first variation of E about u* as 

= / {u* x (j) x — a 2 u*(j) — 4>cosx + X(j)) dx , 
} Jn 



where 

A = — 



^(2E(u*) + J u* cos xdx^j. 



M 

Let be a smooth test function supported in P(u), and set 

:«= [ max 1^11 ) > . 

x&up P (p U*(X) 



By construction, u £ (x) > for |e| < e > and therefore the first variation of E along u e must vanish. 
Since this holds for every smooth test function supported on P(u*), the Euler-Lagrange equation holds 
there. 

Similarly, the first variation is nonnegative for every nonnegative test function </>, because posi- 
tive values of e yield admissible competitors u e on Cm- This means that u* satisfies the variational 
inequality 

u xx + + cos x < A onQ 

in the sense of distributions. Taking = 1, we that A > a 2 |^ > 0, and by considering nonnegative 
test functions supported on Z(u*), we obtain the claimed inequality for A. 

To see that the contact angles are zero, assume that u*{r) = 0. The variational inequality implies 
that u*(x) — ^rr(x — r) 2 is concave in x. In particular, the graph of u* lies below a support line at 
x = r. Solving for u*(x), we see that 

< u*(x) < — - r) 2 + b(x - t) , 

where b is the slope of the support line. It follows that 6 = 0, and u* is differentiable at r with 

u*(t) = 0. □ 

Our next goal is to determine the minimizers of E on Cm explicitly. At first sight, this appears to 
be a simple matter of minimizing over the two free parameters in the general solution of the Euler- 
Lagrange equation. This solution is given by 

u{x) = + u°(x) + A cos(ax) + B sin(ax) , (2.2) 



a 2 



where 



u°(x) = \ 2 ' (2.3) 

[ 13^2 COS X , Q^l. 

A moment's consideration shows that the minimizer cannot be a strictly positive function given by 
Eq. (Eia . unless A = B = and M > (in which case A = ^-). If u vanishes somewhere in 

VL, then A depends implicitly on M through a nonlinear equation. The number of components of the 
positivity set is another unknown, the constants A and B may differ from component to component, 
and the contact points that form the boundary of the positivity set contribute additional free parameters. 

We reduce the number of parameters by observing that minimizers of E on Cm are necessarily 
symmetric decreasing on [— ir, ir] about x = 0. To see this, let w# be the symmetric decreasing 
rearrangement of u 11261 Section 3.3]. By definition, each sub-level set {i G O u#(x) > s} is an open 
interval centered at x = that has the same measure as the corresponding set {x E f2 | u(x) > s}. 
Classical results ensure that u* G Cm, and that 

| \ L 2 = | \u\ |x2 , I \uf\ \ L 2 < | \u x \ \ L 2 , / u*cosxdx> / ucosxdx, (2.4) 

Jn Jn 

see (221 [31. It follows that 



E(u*) < E\ 



u 
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Figure 3: The energy minimizer for a = 0.5, 1, 2 with initial data uq = 3 (left). Mass versus contact point for 
a = 2,1,0.5 (right). 

If u is a minimizer, then E(u#) = E(u), and in particular, the last inequality in (|2.4I) must hold with 
equality. Since the cosine is strictly symmetric decreasing, this forces u to be symmetric decreasing as 
well $M Theorem 3.4]. 

It is now easy to determine the minimizer by solving the Euler-Lagrange equation with zero contact 
angle boundary conditions on a symmetric interval (— r, r). The remaining three parameters are the 
coefficient A, the Lagrange multiplier A, and the contact point r. The next lemma will be used to show 
that the positivity set of a minimizer grows with its mass. 

Lemma 2.3 (One-sided derivatives.) Assume that u* minimizes the energy on Cm- IfP(u*) = (— r, r) 
for some r G (0, n], then u* x (t-) > and A > cosr. 

Proof. Since u* satisfies the Euler-Lagrange equation by Lemma [2T2l it has bounded derivatives of all 
orders on (— r, r). We analyze the sign of the first non- vanishing one-sided derivative of u*. By sym- 
metry, it suffices to consider the right endpoint at x — r. The positivity of u* implies that u* x (tJ) > 0. 

Suppose that u* xx (tJ) = 0, then u xxx (tJ) < 0. On the other hand, by differentiating Eq. (12.11) . we 
obtain u xxx (r_) = sinr > 0. This leaves only the possibility that r = n. Differentiating once more, 
we obtain u* xxx (it-) = —1, which is the wrong sign for u* to have a minimum at it. It follows that 
u* xx (tJ) > 0. The claimed inequality for A follows from Eq. (I2.ll ). □ 

The following theorem summarizes our results, see Figure 3. 

Theorem 2.4 (Description of the energy minimizers.) Let E be the energy functional in Eq. (\1.6\) 

for some a > 0. For each M > 0, E has a unique nonnegative minimizer of mass M. The minimizer is 
strictly symmetric decreasing on its positivity set. It depends continuously on M in C 1,1 , and increases 
pointwise with M in the sense that for any pair of minimizers u\, u\ of mass Mi, M 2 , 

Mi < M 2 =^ u\(x) < u* 2 (x) , x e P(ul) . 

Ifa<l and M(l — a 2 ) > 2tt, then the minimizer is strictly positive and given by 

M 1 

u*(x) = - — h -cos x, x <E Q . (2.5) 

Z7T 1 — a 2 
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Otherwise, it is given by 



u*(x) = u°(x) + Acos(ax) — u°(t) — Acos(ar) , |x| < r (2.6) 

and vanishes elsewhere. Here, r is a smooth increasing function of M with r • maxja, 1} < ir, the 
function u° is given by Eq. (\2.3\) . and A = - • 



Proof. Fix a > and M > 0. By Lemma I2TU there exists a minimizer u* of mass M. If u* is strictly 
positive, then Eq. (12.21) holds for all x E fi. Since positive minimizers are symmetric about x — 0, 
smooth, periodic, and have mass M, we conclude that a < 1 and Eq. (12.51) holds. In order for u* to be 
nonnegative and symmetric decreasing we must have M (1 — a 2 ) > 2n. In that region, u* is clearly 
strictly increasing in M. 

If, on the other hand, the positivity constraint is active, then u* is symmetric decreasing on some 
interval (— r, r) and vanishes for |x| > r. By Lemma [2T2l u* E C 1 ' 1 (fi) and ti*(±r) = 0. On (— r, r), 
it* is given by Eq. (12.21) . Since w* and u° are even, 5 = 0. The Dirichlet condition at r yields 

A = -a 2 (A cos(ar) + u°(r)) , 

the Neumann condition determines A, and we find that Eq. (12.61) holds. We denote this function by 
u*(x; r). If r • maxja, 1} < n, we claim that : r) is indeed nonnegative, symmetric decreasing 
in x, and strictly increasing with r. To see this, we differentiate Eq. (12.61) . and use that u*(r; r) = to 
obtain 

— ■ a sin ar = -< r (r; r) = < x .(r; r) > . 
ar 

By the chain rule, and using once more that u* x (r; r) = 0, we have 

dA 

u*(x; t) = —— ■ (cos(ax) — coslar)) > for |x| < r . 
dr 

Since u* vanishes identically when r = 0, this confirms that it is positive and strictly symmetric 
decreasing for |x| < r. 

The mass of u* is given by M(r) = u*(x; r) dx. We use that u*(t; t) = to compute 



dM dA 



—— [ (cos(ax) — cos(«r)) dx > . 
dr J^ T 



dr dr 

and infer that we can solve for r = r(M) as a strictly increasing smooth function of M. By the chain 
rule and the inverse function theorem, 

d 7 u*(x; r(M)) = ™(^\- > Q . 



dM ' J^ r (cos(ax') — cos(ar)) dx' 

It remains to determine the ranges where Eq. (12.51) and (12.61) hold. For a < 1, the energy minimizer 
on Cm is unique by the strict convexity of E. If M > y^2, the minimizer is given by Eq. (12.51) . 
For smaller values of the mass, we use instead Eq. (12.61 ), and compute that M — > as r — >■ and 
M — >• as r — )• 7r_. Continuous dependence on M follows, since Eq. (12.61) agrees with Eq (12.51) at 

M — 2 " 



1-a 2 
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When a > 1, the positivity constraint is always active, because E is not bounded below without 
it. Therefore, the minimizer is given by Eq. (12.61) on some interval (— r, r). For a = 1, we must have 
r < it, because the particular solution u°(x) = — ^xsinx from Eq. (12.31) cannot be continued as a 
differentiable 2n -periodic function across x = it, in violation of Lemma 12 .2[ It is easy to check from 
Eq. (|2.6I) that Af — > as r — > and Af — >• oo as r — > it. For a > 1, we have that necessarily ar < 1, 
since otherwise the function defined by Eq. (12.61) fails to be symmetric decreasing. Since Af — > as 
r — > and Af — > oo as r — > oT^it, the theorem follows. □ 

3 Steady states 

In this section, we investigate the relationship between steady states of Eq. (11.41) and critical points of 
the energy in Eq. (|1.6I) . For a < 1 and n > 1, we will show that the global minimizer u* determined 
in Theorem 12.41 is the unique point in Cm where the energy dissipation defined in Eq. (|1.7I) vanishes, 
but for a > 1 there are additional steady states. We start with some definitions. 

Definition 3.1 Let u £ Cm such that E{u) < oo, and let P(u) be its positivity set. 

• u G H 2 {Vl) is a strong steady state of Eq. (11.41) if for every smooth 2n -periodic test function (p 

{u xx + a 2 u + cos x) ■ (u n (f) x ) x dx = ; 

• u is an L 2 -critical point of the energy on Cm if every differentiable curve 7 : (— £ , £0) - > Cm in 
L 2 with 7(0) = u satisfies 

£( 7 (e)) > E(u) - o(\\j(e) - u\\ L 2) , as e ->• . 

Strong steady states are time-independent strong solutions in the sense of Eq. (|1.5I) . As elements 
of H 2 , they are of class C 1 ' 1 and meet dry regions with zero contact angles. The reason why we define 
critical points in the L 2 -topology rather than in H 1 is that the boundary of Cm in if 1 , which consist 
of nonnegative functions of mass Af that vanish at some point on f2, contains many curves that are 
differentiable in L 2 but not in H 1 . By analogy with the subdifferential in convex analysis, we ask only 
for a lower bound on the energy difference because E is lower semicontinuous, but not continuous, on 
Cm with the L 2 -norm. 

The following two lemmas relate these notions to the Euler-Lagrange equation. 

Lemma 3.2 (Steady states with zero dissipation.) Let u 6 C M , and define the dissipation D(u) by 
Eq. di.7D . If either 

• u satisfies Eq. (\2.1\) on each component C ofP(u) with a constant A = A(C) and with u = u x = 
on dC, 

or, equivalently, 

• u E Hf oc (P(u)) n ff 2 (fi) with u = u x = 0on dP(u) and D(u) = 0, 
then u is a strong steady state. 
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Proof. We first prove the equivalence of the two conditions. Assume that u G Cm satisfies Eq. (12.11) on 
each component of P(u). Since u xxx = —a 2 u x + sinx G H 1 (P(u)), it follows that u G Hf oc (P(u)) fl 
H 2 (Vt), and D vanishes. Conversely, if D{u) = then u xxx +a 2 u x — sin x vanishes in L 2 oc (P{u)). This 
means that u xx + a 2 u + cos x is locally constant on P(u), i.e., u satisfies Eq. (12.11) on each component 
C of P{u) with a constant A = A(C). 

Let {Cj} be the collection of connected components of P(u). If u solves Eq. (12.11) on each Cj with 
some constant \j and is a smooth 27r-periodic test function, then 



/ (u xx + a 2 u + cosx) • (u n (j) x ) x dx < / Aj • (u n (f) x ) x dx 
Jp(u) „• 1 Jc*, 



0, 



'p(«) 

showing that u is a strong steady state. □ 

Although strong solutions need not be regular enough to justify differentiating the energy, for n > 1 
it is not hard to show (by arguments analogous to [IT71 Lemmas 1 and 2]) that the dissipation vanishes 
in all strong steady states of Eq. (11 .7b . The next lemma generalizes the description of the minimizers 
in Lemma [2T2l 

Lemma 3.3 (Characterization of critical points.) A function u G Cm is an L 2 -critical point of the 
energy if and only if there exists a A G R such that u solves the Euler-Lagrange equation d2.il) with 
this value ofX on every component of P(u) and u = u x = on dP(u). 

Proof. Suppose that u G C M solves Eq. (12.1b on P(u), and that u — u x — on dP(u). For v G C M , 
we compute the directional derivative 



d_ 

ds 



E((l — s)u + sv) = / {u x ■ (v — u) x — (a 2 u + cosx) • (v — u)} dx 



A ■ (v — u) dx — I cos x • (v — u) dx 

P(u) Jz(u) 

v ■ (A — cosx) dx . 

Z(u) 



Since E agrees with its second order Taylor expansion about u, it can be written as 

E(v ) = E(u) + I v ■ (A — cosx) dx + - f {(v x — u x ) 2 — a 2 (v—u) 2 } dx . (3.1) 
Jz(u) 2 J n 

Let 7 : (— e Q , e ) — > Cm be a differentiable curve through u in L 2 . Writing 7(5) = u + s^'{e) + o(e) 
in L 2 , we see that the nonnegativity of 7(e) implies that 7'(e) vanishes almost everywhere on the zero 
setZ(w). By Eq. (13TTT) . 

£(7(e)) - -B(w) > / 7(e) • (A - cosx) dx - - f a 2 (j(e) - u) 2 dx = o(e) , 
J z{u) 2 J n 

showing that u is an L 2 -critical point. 

Conversely, assume that u is an L 2 -critical point. By considering 7(e) = jp^u + e<p), where 4> 
is a smooth function with support in P(u), we see that u satisfies the Euler-Lagrange equation on 
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P(u). Since u G H l , it is continuous and vanishes on dP(u). We need to show that u x also vanishes 
on dP(u). 

Consider a connected component C of P(u), and let £ be its length. By Rolle's theorem, u x 
vanishes somewhere on C, and by the Euler-Lagrange equation, sup^ \u x \ < £ ■ (A + 1 + ||w||oo)- In 
particular, u is Lipschitz continuous on Vt. We claim that u has one-sided derivatives at every point 
r with u{t) = 0. If r is the limit of an increasing sequence of zeroes of u, then it follows from the 
above estimate that u x {rJ) = 0. Otherwise, if r lies on the right boundary of a component C, its left 
derivative exists because u solves the Euler-Lagrange equation on C. Similarly, u x {j + ) exists, and 
vanishes unless r is the left endpoint of a component of P(u). 

Let be a smooth test function on f2, and consider variations of the form 



M 

u £ {x) = — u(x + e<p(x)) 



M £ 



u(x + £(f>(x)) dx . 



Since u G H 1 , the curve 7 : u u £ is differentiable in L 2 with 7'(0) = 0m 2 
chain rule 



We compute with the 



d_ 

de 



E(u e 



£=0 



-ui H — a 2 u 2 + u cos x ) dx — 
2 2 



bu sinx dx — X—M £ 

de 



e=0 



where 



\ = -—[2E(u) + 



I! COS ./• (/./• ) , —M £ 

de 



e=o 



(f)U x . 



Setting the first variation equal to zero yields the (weak) Beltrami identity associated with Eq. (12.11) 
(see, for example lfT51 Theorem 2.8]). We next write P(u) as the union of its connected components 
Cj and integrate the first integral by parts on each component. (The number of components may be 
finite or countable.) Using that u satisfies Eq. (12.11) and vanishes on dP(u), we obtain 



d_ 

de 



E(u £ 



£=0 



dCi 



Since u is an L 2 -critical point, this expression vanishes for every smooth function on fi. By concen- 
trating <fr at one point r G dP(u) , we conclude that | u x (r_ ) | = \u x (r + ) | . 

It remains to show that u x (t+) = 0. Suppose for the contrary that u x {t + ) = —u x {rJ) = a > 0. 
Then r is the common boundary point of two components of P(u). Let / be an open interval that 
contains r but no other zeroes of u, and consider the variation 



u £ (x) 



M 



u(x) , 



x G J, 
otherwise 



M £ 



u(x) dx + / i/j £ (u(x)) dx . 



where e < min xg a/ u(x), and ip £ is defined for y G [0, 00) by 

4> £ (y) = 



y ■ 



y > \ £ \ > 

e > 0, y < e, 
e < 0, y < -e . 



max{2y + e, 0} , 

Then 7 : e 1— > u £ defines a curve in Cm that is differentiable in L 2 with 7(0) = u, and we have 



7 '(0) = 



£ = 



a 

2 ' 
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contradicting the assumption that u is an L 2 -critical point. □ 

Combining Lemma 13.21 with Lemma 13.31 we see that the profile of a steady state where the dis- 
sipation vanishes is the sum of profiles of critical points of the energy whose supports are mutually 
disjoint. We next show that all these profiles are symmetric. 

Lemma 3.4 (Symmetry.) If u is a positive solution of Eq. (\2.1\) on some interval C with boundary 
values u = u x = 0, then it is symmetric about x — 0. 

Proof. Let u be a positive solution of Eq. (|2.1I) on C = (r 1; r 2 ). Then u is given by Eq. (12.21) . For 
a = 1, the boundary conditions u(tj) = u x (rj) = read 

A cos(tj) + B sin(rj) + A — —x sin(rj) = , 

-A sin(r i ) + B cosfo) - - mnfa) - -x cosfa) = . 
After eliminating B, we see that cos r, for j = 1,2 both solve the quadratic equation 

x 2 - 2Xx - c = , (3.2) 

where the constant c is determined by and A and A. For a ^ 1, the boundary conditions can be 
expressed in the single complex equation 

e -*»a (A + iB) = - ^ - t^— 2 cosfa) + - nn(r y ) . 

a 2 1 — a 2 a (I — a 2 ) 

By considering the square modulus A 2 + _B 2 , we see again that cosTj for j = 1, 2 both solve the 
quadratic equation Eq. (13.21) . with a constant that depends on A, B, a, and A. 

By Vieta's theorem, the two roots of Eq. (|3.2I) add up to 2A. Since A > cosTj by Lemma [2721 
we conclude that costx = cosr 2 . This leaves two scenarios: The first is that C = (— r, r) for some 
r G (0, 7r), and u is given by Eq. (12.61) . The other scenario is that C = Q \ [— r, r] for some r G [0, 7r), 
and u is the 27r-periodic continuation of 

u(x) = u°(x) + A cos(a(x + 7r)) — u°(t) — A cos(a(r + n)) , r < x < 2n — t , (3.3) 

with coefficient Afr) = — ^t-V^t- In both cases, u is symmetric about x = 0. □ 

We conclude from Lemmas 13 .2H3 .41 that Eq. (11.41) has four types of strong steady states with zero 
dissipation, see Figure 4. The first three are L 2 -critical points of the energy, but the last generally is not. 
While critical points are typically isolated, the fourth type of steady states can form a one-parameter 
continuum in Cm- 

I. Smooth films are positive except possibly for touchdown zeroes. They occur as energy minimiz- 
ers if a < 1 with M(l — a 2 ) > 2n, and as saddle points if a > 1 with M(a 2 — 1) > 2n. If a is 
not an integer, they are symmetric about x = and given by Eq. (|2.5I) . If a is an integer k > 1 
and M(/c 2 — 1) > 2n, there are also non-symmetric solutions of the form 

, , M 1 _ . , 

u[x) = — cos x + A cos kx + B sin kx 

K J 2tx k 2 -l 

with A and B small enough. 
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Figure 4: Types of steady states with zero dissipation (from left to right): Smooth film, hanging drop, sitting 
drop, and two-droplet steady state. If a = 1, the unique steady state is the energy minimizer, which has the shape 
of a hanging drop. If a < 1, the unique energy-minimizing steady state can take the shape of a hanging drop 
(for small values of the mass) or a smooth film (for larger values of the mass). If a > 1 the energy minimizer is 
always a hanging drop. For larger values of the mass, there are multiple steady states that include sitting drops 
and two-droplet steady states. 

II. Hanging drops are given by Eq. (12.61) on an interval (— r, r) with < r < tt, and have a 
dry region at the top of the cylinder. Among them are the energy minimizers for a < 1 with 
M(l - a 2 ) < 2n, as well as for a > 1. 

The other two types occur only for a > 1: 

III. Sitting drops are given by Eq. (13.31) on an interval (r, 2n — r) with < r < it, and have a dry 
region at the bottom. They are L 2 -critical points, but never minimizers of the energy. 

IV. Two-droplet steady states are the sum of a hanging and a sitting drop whose positivity sets are 
disjoint. They do not correspond to L 2 -critical points of the energy unless the value of A from 
Eq. (12.11) agrees in the two droplets. 

The next theorem is illustrated in Figure 5. 

Theorem 3.5 (Uniqueness of steady states with zero dissipation.) Let M > and a > be given. 
If a < 1, then the global minimizer u* is the unique strong steady state ofEq. A1.4\) on Cm with zero 
dissipation. For a > 1, there exists for each M > an energy level E% > E{u*) such that u* is the 
unique strong steady state with zero dissipation in the sub-level set {u G Cm \ E{u) < Ei}. 

Proof. By Lemmas l3.2H3.4L every strong steady state where the dissipation vanishes is the sum of one 
or two critical points with disjoint positivity sets. For a < 1, the convexity of the energy implies that 
the unique minimizer u* (determined in Theorem 12.41) is the unique critical point. For a > 1, we use 
that the minimizer is unique and given by a hanging drop, and then use Lemma 13.41 to see that it is 
isolated within the set of steady states with zero dissipation. □ 

4 Convergence to equilibrium 

In this section we will prove our main result, that global strong solutions of Eq. (11.41) with a < 1 
converge strongly in H 1 to the unique energy minimizer of the same mass. We have already described 
these steady states in the previous sections. 
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Figure 5: Energy levels of critical points. Dashed line: a = 1. For each given mass, the energy minimizer 
is the unique critical point. Its energy decreases strictly with the mass. Solid line: a = \[2. If the mass 
exceeds a certain threshold, there is in addition to the global minimizer a pair of saddle points. Note that the gap 
between the minimal energy and the energy of the saddles appears to be minimal when the saddles first appear, 
at M = 2ir. Although the energy levels of the two branches cross near M = 8, this is not a bifurcation point. 

As described in the introduction, a global strong solution of Eq. (11.41) is a nonnegative function 
in L 2 OC ((0, oo), H 2 (n)) that satisfies Eq. (11.51) for every smooth test function with compact support in 
(0, oo) x fi. We consider only strong solutions that additionally satisfy a linear bound on the i/ 2 -norm, 



To convince the reader that this class of solutions is not empty, we paraphrase the existence theory of 
Bernis and Friedman for n > 2. We note in passing that the method has been extended to the entire 
range n > 0, and that it implies much stronger existence and regularity results |[4ll7l[T4"l|. 
The basic version of the methods exploits that the entropy 



formally decreases under Eq. (11.31) with n > 2. Note that S(u) < oo implies that u can vanish 
only on a set of measure zero. Strictly speaking, S is not an entropy for Eq. (11.41) , because it may 
increase as well as decrease along solutions, due to the presence of the long-wave instability and the 
gravitational drainage term. These terms can be accommodated with a technique of Chugunova, Pugh 
and Taranets [14J, as follows. Along smooth solutions, 




(4.1) 



with some constants A, B, and the energy inequality 




(4.2) 





(4.3) 
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where c n = (n — 2)(n — 1). Integrating over time and bounding the second integral with Eq. (12.11) . one 
obtains 



S(u(-,T)) + c n I [ u 2 xx dxdt<S(u(-,0)) + KT, (4.4) 
Jo Jn 

where K depends only on the energy and the mass. In particular, classical solutions satisfy Eq. (14.11) . 
A similar computation shows that the energy inequality (14.21) is an identity for classical solutions. 
Strong solutions that satisfy Eqs. (14.11) and (|4.2I) are obtained by regularizing 

u t + d x [f e (u) d x (u xx + a 2 u + cosx)] = , x eQ, 

for example with f e (z) = z n + e |0. The regularized equations are known to have unique smooth 
solutions that exist for all times t > 0. For these solutions, the energy inequality (|4.2I) holds with 
equality, and the entropy inequality (14.31) holds for the corresponding entropy functional with the the 
same constants c n and K. By the usual compactness arguments, there exists a sequence u £j that 
converges uniformly on compact time intervals to a nonnegative strong solution u on (0, oo) x Cl, 
which is smooth wherever it is positive. Moreover, u xx converges to u xx weakly in L 2 (0,T) x Vt 
for each T > 0, and u xxx converges to u xxx weakly in Lf oc (P(u)). In the limit, the energy and 
entropy inequalities in Eqs. (14.21) and (|4.4I) remain intact, because the double integrals are weakly 
lower semicontinuous due to their convexity in the highest derivative. 

We turn to the long-time behaviour of solutions. Lyapunov's principle says that a dissipative dy- 
namical system should converge towards the set of critical points of the energy. Thin film equations 
present two difficulties: Since well-posedness has not been settled, we are not sure how to view them 
as dynamical systems, and lower bounds on the dissipation are not easy to obtain. 

The next lemma bounds the distance from the minimizer in terms of the energy. Since energy 
decreases along solutions of Eq. (|1.4I) . the lemma implies that implies that u* is dynamically stable in 
the sense of Lyapunov. 

Lemma 4.1 (local coercivity.) Let a > and M > be given, and let u* be the energy minimizer on 
Cm obtained in Theorem \2.4\ Then 



sup{d H i(u,u*) ueC M ,E(u) <E(u*) + AE} — ^0 (AE — > 0) . 
For a < 1, we have the explicit estimate 

( 2AE \ 

d H i{u,u*) < [y^) (4 - 5) 

for all u with E(u) < E(u*) + AE. 

Proof. We argue by contradiction. If the conclusion fails, then there exists a minimizing sequence {uj} 
such that inf j (uj, u*) > 0. By Lemma [2TT1 a subsequence converges weakly in H 1 and strongly in 
P? to the minimizer u*. But since the minimizer is unique, and all convergent subsequences have the 
same limit, the entire sequence converges. By the expansion of the energy in Eq. (13.11) . 



d H i(uj,u*) = I 2AE(v,j) — 2 / Uj ■ (X — cosx) dx + a 2 / (uj — u*) 2 dx ] 
V Jz(u*) Jn J 

-> (j -> oo) , 
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contradicting the choice of the sequence. For a < 1, the bound in Eq. (14.51) follows immediately from 
the observation that the linear term in Eq. (|3.1I) is nonnegative, see Lemma [2T2l □ 



For a = 1, one can take advantage of the positivity of the linear term in Eq. (13.11) to obtain an 
explicit estimate of the form dn^(u, u*) < Ci(AE) 1 / 2 + c 2 AE, where the constants C\ and c 2 depend 
on the mass. We next construct a sequence of times along which the dissipation goes to zero. 

Lemma 4.2 (Construction of a weakly convergent sequence.) Let u be a global strong solution of 
Eq. < \1.4\) that satisfies inequalities A4.1\) and A4.2\) . and let E be its energy at time t = 0. There exists 
a sequence of times tj — > oo such that 

sup ||i4 a . a .(-,£,-)|| 2 < oo , lim D(u(-,tj)) = . 

j j^oo 

Proof. Eq. (14.11 ) implies that the set 

C 1 = {te(Q,T)\\\u xx (;t)\\l>4B} 

has measure bounded by /u(Ci) < ^ + j. Similarly, for every e > 0, Eq. (14.21) implies that 

C 2 = {te(0,T)\D(u(.,t))>e} 

has measure bounded by //(C 2 ) < E °~^ u *\ it follows that for T > | + 4 ( g °- g ("*)) > we can find 
t G [j, T] that lies neither in C\ nor in C 2 . The sequence i, is constructed by taking sequences Sj — > 
and Tj oo. □ 

We combine this lemma with the stability result from Lemma |4~T1 to show that u* is in fact asymp- 
totically stable. 

Theorem 4.3 (Asymptotic stability.) Let u be a global strong solution of Eq. of mass M and 
initial energy Eq constructed by the method ofBernis and Friedman, and let u* be the energy minimizer 
on Cm- If ct > 1, assume in addition that the sub=level set {E < E Q } in Cm contains no other steady 
states with zero dissipation. Then 

lim d H i t), u*) = . 

t— >oo 

Proof. Let {tj} be the sequence of times constructed in Lemma |4T2| Since {u(-,tj)} is uniformly 
bounded in H 2 , there is a subsequence (again denoted tj) that converges weakly in H 2 and strongly in 
H 1 to some limit v e Cm- We want to show that v = u*. 

For S > 0, consider the set Ps(v ) = {x E Q \ v(x) > 5}. Since u(-,tj) converges uniformly to v, 
we have that u(x, tj) > | on Ps(v) for j sufficiently large, and it follows that 

/ {u xxx {-,tj) + a 2 u x (-,tj) -sinxfdx < ^D(u(-,tj)) -)■ 0. 
Jp s (v) o 

Since we already know that u x (-,tj) converges to v x strongly in L 2 , this means that u xxx (-,tj) con- 
verges to sinx — a 2 v x strongly in L 2 (Ps(v)). The limit agrees with v xxx , and we see that 



v 



in 



+ a 2 v x — sinx = on P$(v ) . 
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Since 5 was arbitrary and v £ H 2 by construction, it follows from Lemma I3T21 that v is a strong steady 
state of Eq. (11.41) . Since E(v) < E(u(-, 0)) and D{v) = 0, we conclude with Theorem [331 that v = u*. 

We next observe that E(u(-,tj)) — > E{u*) by the continuity of the energy in H 1 . Since the energy 
decreases monotonically along solutions by Eq. (14.21) . we have 

lim E{u{-,t)) = E(u*) , 

t— >oo 

and the claimed convergence follows with Lemma |4~T1 □ 



5 Rate of convergence 

In the final section, we consider the rate of convergence to steady states. We will show that strictly 
positive energy minimizers are exponentially attractive, while steady states that have zeroes can be 
approached at most at a polynomial rate. The reason is that the entropy inequality in Eq. (14.41) limits 
the rate at which the solution can converge to zero on a subset of fi. To obtain the strongest lower 
bound, we will use Kadanoff 's entropy 

S(u) = [ u~ n+ ^dx. (5.1) 
Jn 

One can verify by direct calculation (involving repeated integration by parts flU [20j) that classical 
positive solutions of the thin-film equation Eq. (11.31) with n ^ | satisfy 



_ x dS{u) 



dt 



J u 2 u x u xxx dx = —4 J it 2 ^(it 2 ) D j dx < , (5.2) 



where c n = (n — |)(n — |). This is a special case of Eq. (2.13) in H|. For Eq. (11.41) . Kadanoff's 
entropy can grow at most linearly with time: 

Lemma 5.1 (Entropy bound.) Fix n > |, and let S be given by Eq. {\5.1t Let u £ C M be an initial 
value of finite energy and entropy. Then the global strong solution of Eq. (\1.4\) constructed by the 
method ofBernis and Friedman satisfies 

S(u(;t))<S(u(-,0)) + K t, (5.3) 

where K depends on the mass and the initial energy. 

Proof. If u is a positive classical solution of Eq. (|1.4I) . we differentiate the entropy and integrate by 
parts to obtain 

_,dS(u) f _l , 2 f _1 2 , f 1 

c n — - — = / u 2 u x u xxx ax + a / u ^u x dx — \ u 2 u x sin x ax , 
Jn Jn Jn 

where c n = (n — \){n— \). The first summand we rewrite with the help of Eq. (15.21) . The second 
summand we integrate by parts 

= 2 /«(«»).,««, 
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and combine it with the first by completing the square. This produces a remainder term of the form 
j J^tt? dx. The third summand we integrate by parts as well. We arrive at 



,dS(u) f i / . i. a 2 i\ 2 , a 4 , , , , _ 
c_ — : — = — / u 2 2 M2 u2 dx H — - / «2ai + 2 / «2 cos x ax . 



tft in V 2 ; 4 

By Lemma [2TTI the last two integrals are bounded by a constant that depends only on the mass and the 
energy. Integrating along the solution, we see that Eq. (15.31) holds for classical solutions. 

By the same computation, Eq. (|5.3I) holds for the solutions of a suitably regularized equation with 
the correspondingly regularized entropy and with the same constant K . Since the strong solution is a 
uniform limit of such solutions, the entropy converges as well, and the claim follows. □ 

Theorem 5.2 (Bounds on the rate of convergence.) Consider Eq. ( 17.41 ) with parameters n > and 
a > 0, and set f3 = n — |. Let u be a solution of mass M that satisfies the energy and entropy 
inequalities in ( 14.21) and ( 15.31 ). Assume that u converges in 77 1 to the energy-minimizing steady state 
u* of mass M. 

• If n > | and u* vanishes on a set of positive length L, then 



d H i{u{;t),U*) > 4=- ( 
\/7T V 



\pn \S + K t 



• ifn>2 and u* vanishes quadratically at a point, then there exist positive constants K\ and K2 
(depending on the initial energy and entropy) such that 

d m {u{-,t),u*) > (Kt + Kityifc; 

• if a < 1 andu* is strictly positive, then 

d m (u(;t),u*)<K 3 e-^ 
for some constant K 3 (depending on u), where ji = (1 — a 2 )(min u*) n . 
Proof. If Z(u*) has measure L > 0, we estimate 



S(u(-,t)) > L( sup u(x,t)) ^ > L- \\u(-,t) -u*^ p 



Since u and u* have the same mass, we have — < d H i(u(-,t),u*), and the first 

claim follows from the bound on the entropy in Eq. (|5.3I) . 

If u* has a zero of order 7 > i, we consider the interval I of length L centered at that point and 
obtain with the same calculation as for the first case that 

I \u(-, t) — u*\ > sup u(x, t) — sup u*(x) > ( — — - J — 0(1?) 

xei xei \So + Kot/ 

as L — > 0. Choosing L = e ■ (S + K t)~ "t- 1 , we see that for e > sufficiently small 

|K-,*)-u*|| £O0 > (ATi + 7T 2 t)-^ 
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with some constants K 1: K 2 > 0. Setting 7 = 2 and adjusting the constants we obtain the second 
claim. 

If u* is strictly positive, then a < 1 by Theorem l2.4l By Eq. (13.11) . 

E(u{-,t)) - E{u*) = ~ / {{u x - <) 2 - a 2 {u - u*)f dx 

= 7T - q2 )\Hp) ~ U*(p)\ 2 . 

pGZ\0 

Since u converges to u* in H 1 , there exists a time to such that min u(-,t) > for all t > to- At all later 
times, u is a strictly positive, classical solution that can be differentiated as often as necessary. Since 
u* solves the Euler-Lagrange equation (12.11 ), the dissipation satisfies 

— E(u(-,t)) = - [ u n ((u- u*) xxx + a 2 {u - u*) x f dx 
dt J n 

< — (min-u) n - / [d x ((u — u*) xx + a 2 {u — u*))] 2 dx 

Jn 

= — (minu) n -7r p 2 (p 2 — a 2 ) 2 \u(p) — u*(p)\ 2 

pez\{o} 

/minw \ ra . . . NX , ,., 

< - i^-J -2/i (E(u(-,t ))-E(u* )). 

\mmu / 

In the last two steps, we have used Parseval's identity to rewrite the integral in terms of the Fourier 
coefficients of u — u*, and estimated the Fourier multipliers by 

p 2 (p 2 -a 2 ) 2 >(l-a 2 )(p 2 -a 2 ), (p + 0) . 

Since u(-, t) converges uniformly to u* as t — > 00 by Theorem 14.3 [ it follows from Gronwall's lemma 
that E(u(", £)) — E(u*) < Ke~ 2 ^ for some constant K. By Eq. (14.51) of Lemma |4~T1 this implies the 
claimed exponential convergence of g?#i (u(-,t),u*). □ 

If n = I and w* vanishes on a set of positive length, we obtain as in the proof of the first case of the 
theorem yields an exponential bound of the form d H i («(•, t),u*) > Ke^^, where K and /x depend on 
the mass, energy, and entropy of the solution. 

Summary (Pukhnachev's model on a stationary cylinder.) Let ube a global strong solution of 

u t + d x [u 3 d x (u xx + u + cos a;)] = , x G R/(2irZ) 

constructed by the method of Bernis and Friedman, and let u* be the unique nonnegative energy 
minimizer of the same mass M. Then 

lim d H i(u(-, t),u*) = . 

t— >oo 

The minimizer is a droplet with zero contact angles and profile 

u*(x) = — -{x sinx — r sin r) + -(1 — r cot r)(cosr — cos x) , \x\ < r . 
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Figure 6: Evolution of a solution with a = 1, n = 3 and initial data u$ = 1. Time shots of the numerical 
solution at t = 0, 10~ 2 , 10 -1 , 1, 10, 10 2 , 10 3 (left). L 2 -distance of the solution from the energy minimizer. The 
dashed line shows the lower bound from Eq. ( 15.41 ) (right). 



The contact point r is a continuous, strictly increasing function of the mass, with t — > as M — > 
and t — > 7i as M — > oo. If, additionally, 



/ (u(x,t)) 2 dx < S + K t, 
Jn 



then 

dHM , t) ,u-) >-L.(|^)». (5.4) 
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